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Abstract: We present a fit to the 2012 LHC Higgs data in different supersymmetric frameworks 
using naturalness as a guiding principle. We consider the MSSM and its D-term and F-term extensions 
that can raise the tree-level Higgs mass. When adding an extra chiral superficld to the MSSM, three 
parameters are needed determine the tree-level couplings of the lightest Higgs. Two more parameters 
cover the most relevant loop corrections, that affect the /177 and hgg vertexes. Motivated by this 
consideration, we present the results of a five parameters fit encompassing a vast class of complete 
supersymmetric theories. We find meaningful bounds on singlet mixing and on the mass of the 
pseudoscalar Higgs as a function of tan/3 in the MSSM. We show that in the (rattan/?) plane, 
Higgs couplings measurements are probing areas of parameter space currently inaccessible to direct 
searches. We also consider separately the two cases in which only loop effects or only tree-level 
effects are sizeable. In the former case we study in detail stops' and charginos' contributions to Higgs 
couplings, while in the latter we show that the data point to the decoupling limit of the Higgs sector. 
In a particular realization of the decoupling limit, with an approximate PQ symmetry, we obtain 
constraints on the heavy scalar Higgs mass in a general type-II Two Higgs Doublet Model. 
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1 Introduction 

With the 2012 run reaching its conclusion our knowledge of Higgs couplings is considerably improving. 
Some purists might argue that the new resonance [1, 2] is not yet experimentally proven to be the 
Higgs boson, but all evidence goes in that direction and even more so after the recent updates by the 
LHC experiments [3, 4]. In the previous round of LHC results, namely 4.7 fb _1 at y/s = 7 TeV and 5.9 
ftV 1 at ^ = 8 TeV, the errors in the Higgs' couplings were dominated by statistical uncertainties that 
are expected to scale-down with increased luminosity. Nonetheless, these measurements were already 
precise enough to make meaningful statements about broad classes of models [5-9] . After the addition 
of approximately 7 fb _1 at ^fs = 8 TeV, we find it interesting to explore the consequences of Higgs 
measurements in the predictive and theoretically motivated framework of natural supersymmetry. 
Several groups have already fitted Higgs' rates in different contexts, focusing mainly on simplified 
settings in which all the couplings were determined by one or two free parameters [10-14] , concentrating 
also on the cases of light stops or the type-II two Higgs doublet model (2HDM) [5, 10]. Here we 
adopt a different perspective, discussing the implications of Higgs measurements in complete natural 
supersymmctric theories. 

In Section 2 we briefly review our definition of natural supersymmetry. In this context, it is 
necessary to extend the Higgs sector of the MSSM to obtain the observed Higgs mass rrih ps 126 GeV 
without incurring excessive fine-tuning. Therefore we consider D-term and _F-term models that modify 
the 2HDM structure of the MSSM potential. 

We study the new tree-level effects, distinguishing between the cases in which an approximately 
type-II 2HDM structure is preserved and those in which the mixing with extra light states in the 
Higgs sector is sizeable. We find that in general three couplings are enough to parameterize all the 
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tree-level deviations in Higgs rates. Loop corrections to hgg and hyy vertices can also be sizeable, 
since naturalness considerations require stops and Higgsinos to be fairly light. Two extra parameters 
are needed to include these effects. Therefore we perform a five dimensional fit to the data. Although 
naturalness has driven us to consider a certain spectrum and a certain subset of supersymmetric effects 
on Higgs couplings, the fine-tuning constraints are not imposed in fits. The general five dimensional 
framework parameterizes a vast class of complete theories, that may or may not be fine-tuned. We also 
present fit results in several simplified settings, considering the cases in which: (1) only the tree-level 
mixing with the second Higgs doublet plays a role, (2) only tree-level effects are important, but both 
the second Higgs doublet and an extra singlet are responsible for deviations in the couplings and (3) 
only MSSM loop effects from stops and charginos are sizeable. 

The fitting procedure is described in Section 3 together with the data taken as input. In Section 4 
we present the results of the fit and discuss their relevance in light of direct searches performed at the 
LHC. 



2 The Natural SUSY framework 

Deviations in the Higgs couplings are probed by the measurements of the Higgs rates, which are 
determined by partial widths. We denote the modifications to the partial widths as 



i H 2 = r ( fe -^) (21) 



with i = t, V, G, 7, b, t standing for top, massive vector gauge boson, gluon, photon, bottom and tau, 
respectively. In general it is hard to write the r^'s in terms of few supersymmetry breaking parameters. 
However if we restrict to natural theories it is possible to identify a limited number of relevant effects, 
as was shown in [15]. As stated above, we do not impose fine-tuning bounds in the fits, but we still find 
naturalness as a good guiding principle to select models and simplified scenarios in which fit results 
can be interpreted. 

For concreteness, consider at most a A -1 = 10% tuning, where 



2Sm% 



.- I ^ „.2 , 2 , I a |2^ - A 



^J st0 ps = -^Vt [m 2 Q 3 + r< 3 + \A t \ J J log — . (2.2) 



This restricts, in the MSSM, an unmixed stop to be lighter than about 600 GeV. Taking into account 
also the relation between the Higgs mass and the Higgsino mass: —m\/2 = |^| 2 + .... limits one 
chargino to be lighter than 300 GeV. More precisely, the upper bounds obtained from fine-tuning take 
the form [16] 



smj3 / 3 m h 10% 



< + mg, < 8 80 G.V— ^ ^ bg(A/TeV) 125 » eV V ^? ■ (2-3) 



where a t = A t / \J m \ i + m \ an d A is the scale at which supersymmetry breaking effects are mediated 
to the MSSM. At tree level, the bound on /j, translates into a bound on the lightest chargino 



m x ± < y/lMl 2 + < 300 GeV. (2.4) 
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Other Standard Model (SM) superpartners can be much heavier, with masses ranging from a few TeVs 
to tens of TeVs. 

However, raising the Higgs mass to 126 GeV in the MSSM through stop radiative corrections, 
implies a tuning much worse than 10% [17, 18]. To accommodate a 126 GeV Higgs boson, natural 
supersymmetry must be non-minimal, allowing for additional interactions that modify the Higgs sector. 
The new interactions deform the Higgs quartic potential at tree level and change the Higgs couplings 
beyond the significant deviations (relative to the SM) that can already occur in the MSSM. In the 
following we will consider both extensions that leave the type-II 2HDM structure of the potential 
unaltered and shift quartic couplings relative to the MSSM, and models in which the 2HDM relations 
between light Higgs couplings are relaxed. 

Nonetheless, naturalness not only complicates the problem, as seen above, it also limits the size 
of At and /i. This greatly simplifies the problem of determining the r*j's as a function of soft super- 
symmetry breaking parameters. First and foremost, the type-II 2HDM structure of the MSSM Higgs 
potential is not strongly affected by non-holomorphic loop corrections in a natural theory, regardless 
of our definition of tuning. As an example we can take the definition at the beginning of this section 
and ask for a 10% tuning, then the corrections to the bottom Yukawa for tan/3 < 40 (and r& = 0(1) 
at tree-level) are negligible [15]. So in the natural MSSM and its D-term extensions, Higgs couplings 
can be expressed in terms of only two parameters at tree-level, as in a general type-II 2HDM. Adding 
a singlet chiral superfield, as for example in the NMSSM, increases the number of parameters to three, 
if we assume all couplings to be CP conserving. The additional degree-of-freedom parameterizes the 
mixing between the new CP even state and the lightest scalar Higgs: = (S\h). 

We have seen that naturalness also prefers stops and charginos to be fairly light, and for this reason 
we consider their loop contributions to the Higgs to digluon/diphoton partial widths, introducing two 
new parameters. They can be chosen as Sr G , characterizing the stop contribution to the dimension- 
5 Higgs-gluon-gluon coupling and 5r* , characterizing the chargino contribution to the dimension-5 
Higgs-photon-photon-coupling. It was suggested that light staus can be responsible for the enhance- 
ment of the h — > 77 rate observed at the LHC [19-21], but we do not consider this possibility explicitly 
in what follows. We also neglect sbottom loop corrections and charged Higgs loop corrections. In all 
three cases we expect contributions larger than a few % only in extreme corners of the parameter 
space that imply a tuning much worse than 10% in our definition [15]. In this regions charge breaking 
minima can also be generated [21]. Nonetheless, in the 5 dimensional fits, tq and r 7 are treated as 
independent parameters. Thus we capture any loop effect on these vertexes. 

In terms of the five parameters discussed above, the partial width modifiers are given by 

r T = r b 

rv = (1 - 4 + r b r t )/ (r t + r b ) ^ 

r G = r t (l + 5r l G ) 

r 7 « 1.27r v -0.27r G + Srf , 

where we have chosen and r t to parametrize the 2HDM tree-level effects. Note that r r = rb 
assumes that non-holomorphic corrections to the Higgs interactions with down quarks and leptons 
are small relative to the tree-level value, as discussed above. Here, and in what follows, we always 
neglect Higgs decays to non-SM particles, in particular to neutralinos. In Section 4.1 we review all 
tree-level predictions in more detail, starting with a generic 2HDM with = 0, and later introducing 
a superpotential for the singlet. In Sections 4.4 and 4.5 we study the bounds on loop corrections from 
naturalness and direct searches and compare them with the results of the fit. 
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3 Data and fitting procedure 



In order to constrain the supersymmetric parameters, we take into account all the available Higgs 
channel rates [3, 4]. The rates depend on the product of the overall production cross-section and 
branching ratio for the particular channel. The results are typically reported as a confidence interval 
on the event rate relative to the SM prediction, denoted by (i. We take all the /i's at the two LHC best 
fit values for the mass: m h = 125.8 GeV for CMS and m h = 126.0 GeV for ATLAS. For the Tevatron 
h — > bb rate we take the value of /2 at = 126.0 GeV. In Table 1 we list the analyses relevant for 
our fit, regrouped by channel, and the corresponding values of the /i's. 

In any given model, the signal strengths are determined by the r^'s. For the LHC, the four relevant 
production modes and their respective theoretical dependance are 

• Gluon fusion (gg — > h): a GF /a^ = |r G | 2 , 

• Vector boson fusion (qq — > hqq): <Jvbf/&vbf — \ r v\ 2 , 

• Vector boson associated production (qq — > hV): <tvh/^^h = \ r v\ 2 , 

• Top associated production (gg — > htt): cr ht i/a^ = \r t \ 2 . 

In some cases, a channel can include events from several production modes; for instance, the dijet 
tagged signature in Table 1 is dominantly produced via vector boson fusion, but contains a non- 
negligible contribution from gluon fusion. The collaborations have made public all the numbers needed 
to assess the composition in terms of the physical production modes in the references listed in Table 1. 
In the fit we use the numbers provided channel by channel, but we find a good uniformity between 
different final states and the two experiments, giving roughly 75 — 80% of VBF and 20 — 25% of GF 
in the dijet tagged categories. The relative fraction of the untagged mode are close to be the relative 
cross-section fractions of the different channels. Therefore, where efficiencies are not publicly available 
(e.g. ATLAS untagged WW), we use the ratio between cross-sections to determine the composition 
of the sample. The VH and UH categories can be taken as pure. 
The fits to the data are performed minimizing the % 2 



where Hi(rj) are predicted rates and cr^'s are two-sided errors (depending on the sign of \Xi — jlj). By 
performing the \ 2 test we are assuming that the likelihood functions for fi follow an approximate two- 
sided gaussian distribution and the correlations can be neglected, both of which have been shown to 
be valid approximations [5, 13]. Furthermore we have compared our results with the two dimensional 
fits released by the collaborations, obtaining a good agreement. 

In the following we refer to "preferred" and "allowed" regions. Preferred regions are obtained 
by varying an A-parameters subset of the r^s, while fixing the other parameters to their SM- 
values. Strictly speaking, allowed regions are found including all the r^'s, the possibility of flavor 
non- universality (for instance r\, ^ r T ), vertex structures different from the SM and the presence of 
an invisible width. At this stage, with limited precision in the measurements, we find a theoretically 
inspired (see Sections 2 and 4) five dimensional fit a reasonable approximation of the most generic 
setting. In this framework we obtain lower dimensional confidence intervals by treating the remaining 
independent parameters as nuisances (given a lower dimensional point the other parameters are varied 
to give the best possible fit). 




i— channels 



(3.1) 
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Table 1. Rates relative to their SM value with the respective 68% confidence interval, as measured by the 
ATLAS (top), CMS (middle) and CDF+DO (bottom) experimental collaborations. 



ATLAS 


untagged 


dijet tagged 


VH 


ttH 


h — > 77 


1.80 ±0.50 [22] 


2.7 ± 1.3 [22] 


X 


X 


WW 


l-34i° ; 5 5 9 [23] 


X 


X 


X 


ZZ 


l-20±g;i° [24] 


X 


X 


X 


h — y tt 


(*)x 


0.69i°:£[25] 


X 


X 


h^bb 


X 


X 


0.7 ± 1.1 [26] 


X 


CMS 


untagged 


dijet tagged 


VH 


ttH 


h — > 77 


1.48l° ;g[27] 


2-17±^ 3 [27] 


X 


X 


WW 


0.79l° ; 2 2 6 7 [28] 


0.00l°;™ [28] 


-0-211:5 [28] 


X 


ZZ 


0-82±g;|* [29] 


X 


X 


X 


h — > TT 


o.8+l : | [30] 


0.82±° f 5 [30] 


i.il?:? Pi] 


X 


h^bb 


X 


X 


1.33lg:» [32] 


-0.8 + | 2 [33] 


CDF+DO 


untagged 


dijet tagged 


VH 


ttH 


h^bb 


X 


X 


1 fi1+ 074 


X 



(*) We have not included the non-VBF tt categories for ATLAS. The information on the cross 
correlation is not provided making difficult to asses the contamination of the combined result. To 
extract the VBF only rate from the CL S plot we follow the procedure described in [14]. 

4 Predictions and Results 

In this Section we discuss in more detail the predictions of natural supersymmetry and compare them 
with the results of the fit to Higgs data. We consider four simplified scenarios: a type-II 2HDM , a 
type-II 2HDM with a SM singlet, light stops and light charginos both in the Higgs decoupling limit. 
At the end of the section we discuss the more general five parameters fit. 

4.1 The type-II Two Higgs Doublet Model 

If the superpartners are relatively heavy, i.e. > 500 GeV and m x ± > 200 GeV, loops corrections 
to Higgs couplings are at most of the order of 10% in a broad region of parameter space, provided that 
the stop mixing is limited by naturalness [15]. Therefore it is easy to decouple these effects without 
incurring excessive fine-tuning and without the need to single out narrow corners of parameter space. 

With this in mind, we begin by exploring the tree-level corrections in the MSSM and its extensions, 
using the generic type-II 2HDM scalar potential and Yukawa couplings [34, 35], 

-£ = m 2 \H d \ 2 +m 2 \H u \ 2 

+ yl^l 4 + y \H u f + \ 3 \H d \ 2 \H u \ 2 + \4HlH u \ 2 

+ ^ (H d ■ H u ) 2 + (H d ■ H u ) (m 2 2 + X 6 \H d \ 2 + X 7 \H U \ 2 ) + h.c. 
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+ Y t H u ■ Q 3 t c R + Y b H d - Q 3 b c R + Y T H d - L 3 r R + h.c. , (4.1) 

where A-B = e lJ A i B ] , and we adopt the usual definition H u = (if+, if°) T and H d = (H%, H d ~) T '. We 
also take all couplings in the Higgs potential to be CP-conserving. In this setting the Higgs couplings 
to the SM particles depend on two parameters, the ratio of Higgs vevs [34] 

tan/3 =™ (4.2) 

and the mixing angle of the CP-even Higgs states, defined by 

(if \ i-( cos a sin a \ ( Rc if 9 \ 
) =V2[ J . (4.3) 

/i J y — sin a cos a / \ Re if„ / 

where /i is the observed 126 GeV state. At tree-level, r& = r T , 

sin a cos a . . . . . 

n, = 3, r t = -— rv = sin (/3 - a) , (4.4) 

cos p sin p 

and in the MSSM 

sin 2a m 2 H + m\ tan 2a m 2 A +m 2 z . , 

sin 2/3 m H ~ m h' tan 2/3 — m| ' 

Without loss of generality, we can take tan/3 > [36], with Eq. (4.5) implying sin a < 0. We require 
that the top Yukawa does not blow up above the electroweak scale, which imposes the lower bound 
tan/3 > 1 [37]. As discussed in Section 2, demanding loop corrections to the type-II 2HDM structure 
to be negligible provides the additional constraint tan /3 < 40 (that can be relaxed or made stronger 
depending on the level of fine-tuning that we allow in the theory). 

Since the measured rates broadly agree with the SM, the results of the fit, which are shown in 
Figure 1, point to the decoupling limit of the model, £ = a — f3 + ir/2 ps 0. This translates into 
the bound |£| < 0.1 at 95% C.L. Note that the discontinuity in the % 2 is physical and can be easily 
understood by expanding Higgs couplings for small £ 

r b = l-£tan/3 + 0(e 2 ), n = 1 + f / tan/3 + 0(£ 2 ), r v = 1 - ^ + 0(£ 4 ) . (4.6) 

Even for extremely small values of £ the correction to r b can be significant for large tan /3 and depends 
on the sign of £. The one dimensional x 2 wa s obtained by treating tan /3 as a nuisance, so the profiling 
is selecting a tan /3 that can contribute an observable deviation on the side where r b is depleted, while 
for £ < a tan/3 giving the smallest possible deviation is singled out. For this reason the offset between 
the two sides is the same as the offset between the deepest minimum and the SM. 

The result, albeit disappointing from the perspective of finding new physics, allows us to consider 
a particular realization of the decoupling limit in which more interesting statements can be made on 
the mass of the second Higgs doublet. In general it is difficult to translate the fit results for a and 
/3 into a direct mass exclusion or into an exclusion on the Lagrangian parameters. However, if the 
2HDM potential respects an approximate PQ symmetry, or a Z2: H d — > —H d , H u — > if„, advances 
can be made. 

Assuming that the second Higgs doublet is parametrically heavier than fi„, we can expand the 
light Higgs couplings in the Z2 breaking spurion B/M 2 , where 1 

M 2 _ 2 1 A 35 {Hlf _ 2 A 7 (if°) 

M 1 - m 1 H , B - m 12 H . (4.7) 



1 In the following we loosely refer to this expansion as arising from an approximate PQ symmetry, but A5 can be 
large without spoiling our approximation and preserving only the Z2: — ¥ —Hd, H u —> H u . 
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sin a ^=a—/3+7:/2 



Figure 1. Left: x 2 contours corresponding to the 95% confidence level in the (sin a, tan 0) plane (for a type-II 
2HDM). Note that in the MSSM sin a < 0. The black dashed line corresponds to the decoupling limit of the 
MSSM a = P - 7r/2. Right: A* 2 vs. f (= a - /3 + tt/2). The red lines mark the 68% and 95% confidence 
levels heights of the Ax 2 - In black the offset between the SM value of the % 2 an d the 2HDM value at the 
minimum. The discontinuity in the \ 2 is physical and it is illustrated by Eq. (4.6). Note that in the tree-level 
MSSM e > 0. 



and A 35 = A 3 + A 5 . Up to leading order (B/Mi) = 1/ tan/3, thus the expansion converges rather fast 
in 1/ tan/3, e.g., only a O(10%) error for tan/3 = 3. This setting is extremely predictive, giving (in 
the MSSM and many of its extensions where A7 = 0) [38] 



r t = r v = 1 + 0(1/ tan 2 /3) . 



tan 2 PJ J (4.8) 



In this limit, the corrections to rj, depend on a single combination of quartics, A35, and the mass of 
the heavy scalar Higgs mn- Note that if we assume the approximate PQ and some degree of custodial 
SU(2) symmetry, the differences in the masses of the heavy Higgses are small [39]: tuh ~ wia — ttih+- 
Thus, to first approximation mn can be considered the mass of the full heavy doublet. 

The result of the fit is depicted in Figure 2. In the MSSM, A 3 = -{g 2 + (s') 2 )/ 4 « -0.14 and 
A5 = at tree-level, while the loop corrections are small. This leads to a preferred region at 95% C.L. 
with mn > 370 GeV. However the tan /3 expansion clearly introduces an error that in the MSSM goes 
in the direction of relaxing our bound 2 

4m| 12ml 1 „ /mi 1 



! + ^# ^ ITT) ® ( '~~~T 1 — ) ■ ( 4 - 9 ) 

" J tan 2 /3 tan' ;/ 



2 The approximation in (4.9) is valid only at tree-level and the leading term can be obtained from Eq. (4.8) ignoring 
the loop corrections to A2: — > mz (1 + 1/ tan 2 /3). 
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Figure 2. Left: x 2 contours corresponding to the 68% and 95% confidence levels in the (A35,m#) plane, 
where r-,, r t and rv are given by 4.8. The blue dashed line corresponds to the best fit points, the solid purple 
line to the tree-level value of A35 in the MSSM. The band in yellow covers to the possible values of A35 in pure 
D-term models (defined in section 4). 

In the following section we obtain a limit on tua as a function of tan j3 using the full tree-level expression 
in the MSSM together with the most important loop corrections. The effect of taking into account 
the 0(1/ tan 2 (3) variations of the other couplings is much smaller and will be commented upon while 
repeating this exercise in the five dimensional fit case and in the MSSM with the full tan f3 dependence. 

This expansion is also relevant to a vast class of models in which the supersymmetric Higgs sector 
is not necessarily minimal, but it is still well described by a type-II 2HDM at low energy. This is the 
case for D-term models raising the Higgs mass via extra gauge interactions and also for some areas of 
the parameter space of D-term models. 

We define D-term models as the class of theories in which the two Higgs doublets are charged under 
new gauge interactions. We will restrict to gauge theories that do not forbid the /x term, otherwise 
additional D-terms should be present to generate it. In D-term models the bb rate is generically 
enhanced. This descends from the fact that if the charge assignment allows the presence of a /i term, 
the new D-term contributes to the relevant part of the scalar potential 

(|^f -Kf +^) 2 , (4.10) 

with some positive definite coefficients in front. £d contains both gauge symmetry breaking terms 
that can be absorbed by the soft masses m\ and m| and new fields charged under the gauge groups. 
Thus the new effect will be a contribution to A3 that is always negative, resulting in stronger bounds 
on rriH- This is illustrated in Figure 2, by the yellow band corresponding to variation of the D-term 
contribution to the Higgs mass from to — mz- Note that to obtain an observable effect from the 
extra gauge sector, it is necessary to have a large supersymmetry breaking mass for the new gauge 
bosons. If the heavy gauge bosons are integrated out supersymmetrically, there is a shift in the quartic 
<5A 3 ~ [J?/M%r I 40 ]: where fi is limited by naturalness while M v > 3 TeV due to EWPTs [41]. 
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For F-term models the applicability of the approximation in Eq. (4.8) is limited, but we can 
consider a superpotential of the form 



W = \SH u H d + f(S) , (4.11) 

where S is either a singlet under the SM gauge groups or an SU(2) triplet with Y = 0. Integrating out 
S and neglecting v 2 /m 2 s effects (where m s is a supersymmetry breaking mass), the leading observable 
modification to the Higgs potential is 

A MSSM ^ A MSSM + | A |2 ^ (4J2) 

which again shifts A3 by a definite sign, but in the opposite direction with respect to Z?-term models. 
Note that in this case small tan /3 may be preferred to lift the tree-level Higgs mass, 

m 2 h < m| (cos 2 2/3 + 2|> f ' sin 2 2/3 ) (4.13) 
V 9 + (9) J 

and the expansion given in Eq. (4.8) may not be theoretically motivated. Furthermore to leading 
order, in gauge and gaugino mediation, the singlet supersymmetry breaking mass m s in the NMSSM 
vanishes at the mediation scale [42], so it is necessary to introduce some deformations of the simplest 
scenarios to lift the mass above the other soft terms at the electroweak scale. There are nonetheless 
large areas of parameter space, in which the expansion is reliable, that give the correct Higgs mass. 
For instance, including stop corrections with w mj ~ 380 GeV and X t = 0, we get m^ = 126 
GeV for all values of A between 1 and 2, varying tan/3 between 3 and 8 [15]. Considering this range 
for A would correspond to taking A35 between 1 and 4, thus giving a favored region, that for A = 1, is 
already m H > 460 GeV at 95% C.L. 

4.2 Heavy Higgses and m h = 126 GeV in the MSSM 

As discussed above only two parameters are needed in the MSSM, at tree-level, to specify all Higgs 
couplings. Therefore we can express the fit constraints in the (toa, tan/3) plane and compare our 
results with direct searches. To do so, the requirement of naturalness needs to be abandoned to allow 
for the stops' radiative corrections to raise the Higgs mass to 126 GeV. Considering only the leading 
loop correction (which is the only one that survives in the limit fi — > 0) only the H u — H u element 
of the CP even mass matrix receives a correction. The tree-level relations in the previous section are 
correspondingly shifted to [43] 

tan 2a _ m 2 A +m% 



tan2/3 m 2 A - m %+s m y c OS 2p > (4*14) 
m 2 H = —m\ + m\ + m 2 z + Sm 2 , 

where dm 2 can be approximated by the leading stop correction 



x 2 _ 3m t 
dm 



2n 2 v 2 sin 



, M? X 2 / X 2 

log -4 + ^ 1 



2M? \ 6M? 



M{= -{m k +m i2 ) , (4.15) 



but in the following we always fix Sm 2 to the to^ and tan /3 dependent value giving the correct Higgs 
mass 

p 9. m^(m 2 A — m\ + m 2 7 ) - m 2 A m 2 7 cos 2 213 ,.*„s 




Figure 3. Comparison between the fit exclusion in the (rriA,tan/3) plane and the direct exclusions from CMS 
[44] and LEP [45]. The excluded region was obtained profiling the full five dimensional \ 2 ■ Note that the 
experimental collaborations use the m™ ax scenario [46] to set their limits, without imposing rrih ~ 126 GeV. 
Varying the soft supersymmetry breaking parameters might lead to modifications of the observed bounds [47] . 



Before describing the results of the fit we review current LHC searches for MSSM Higgses. Direct 
searches for the neutral MSSM Higgses are currently performed in the cf> — > t + t~ [44, 48], bb(j> — > 
bbfifj, [49] and bbcj) — > bbbb [50] channels. The strongest bound is set by the recently updated CMS 
measurement [44] and ranges from itih > 250 GeV for tan/3 = 5 to run > 700 GeV for tan/3 = 40. 
The sensitivity vanishes below tan/3 = 5. These searches benefit from large tan/3, both thanks to 
new production mechanisms that become important (6 and bb associated production for instance) 
and from the increase in the branching ratio to t + t~. We also include the LEP bound [45] in our 
comparison with Higgs rates, but we do not consider the implications of searches for the charged Higgs 
that are currently not as sensitive as the </> — > t + t~ one plus the LEP constraint. The experimental 
collaborations use the m™ ax scenario [46] to set their limits, and varying the soft supersymmetry 
breaking parameters might lead to modifications of the observed bounds [47] . The effects are stronger 
for large tan /3 mainly due to loop corrections to and can not produce any significant gain in 
sensitivity in the region 3 < tan/3 < 12, where production cross sections become too small and our 
analysis starts to be competitive. Additionally, the collaborations do not impose mt ~ 126 GeV. 
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Ideally, the bounds from [44] should be reinterpreted as a bound on a(pp — > <j})Br{(j> — > t+t~), but it 
is beyond the scope of this work. 

Figure 3 shows the overlay of these bounds with the results of the fit to the Higgs data. In the 
slice 2 < tan/3 < 5 Higgs rates are probing regions of the parameter space not directly accessible to 
CMS and LEP. The fit results in Figure 3 are an anticipation of the five dimensional fit discussed in 
section 4.6 since the excluded region was obtained profiling the five dimensional x 2 j treating r 7 and tq 
as nuisances. This is reflected in the fact that for large tan/3 there is a lower bound tua ^ 250 GeV, 
in agreement with the result of the PQ expansion in Figure 8. 

4.3 Singlet Mixing 

Finally, we consider the effects of the singlet mixing with the two Higgs doublets. We take the 
superpotential (4.11) and do not specify the detailed form of the soft supersymmetry breaking terms 
or extra superpotential interactions. For simplicity we take all new couplings to be CP conserving 
and assume the SM-like Higgs to be CP even and the lightest Higgs state in the theory. With these 
assumptions we do not need to consider exotic decays that could be triggered by approximate Z7(l)'s, 
including the U(1)r, generating light pseudo-goldstone bosons in the Higgs sector. This considerably 
simplifies the problem, leaving the relevant parameters to be contained in the 3x3 mixing matrix for 
the scalar CP even states, 



(h 2 \ 








(s\ 








Sa^x CoiC<pSx 




H d 


(4.17) 


V h J 
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where c x — cosx,s x = smx. In this notation, the tree-level Higgs couplings depend on («,</>,/?), and 
Eq. (4.18) is modified to 

sin cos ex. 

rb = —sin(j) -, r{ = sm<^> -, ry = sin <j> sin (/? — a) , (4-18) 

cos p sin p 

where cos0 = (S\h) measures the amount of singlet in the lightest Higgs. The relations (4.5) are 
only valid in the MSSM, so that tua is no longer determined solely (at tree level) from the mixing 
parameters a and /3. Thus, we do not impose any constraint on a and 0, but still restrict to tan/3 > 1. 
In Fig. 4, we plot A% 2 vs cos 2 </>. Since the presence of singlet mixing will only decrease the observed 
rates, there is no improvement in the fits over the SM and we obtain an upper bound | cos</>| < 0.7 at 
95% C.L. It is worth pointing out that even though this bound was derived assuming only tree-level 
effects are present, it is similar to the bound obtained in the full 5D space, as will be discussed in 
Section 4.6. 

4.4 Stops 

LHC searches, together with naturalness, have stimulated a vast offspring of models of mediation that 
produce light stops, decoupling all other squarks [51-56]. In addition, bottom-up approaches, driven 
mainly by FCNC constraints, still allow stops to be fairly light [52, 57, 58]. It is then worth exploring 
the stop plane in view of Higgs data. 

Using the Higgs low energy theorem [59-61] it is straightforward to obtain the change in the gluon 
fusion rate from integrating out the stops (neglecting Z3-terms) 
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Figure 4. A^ 2 as a function of the singlet mixing parameter cos 2 <j>. The 68% and 95% confidence levels are 
marked in red. In black we show the offset between the \ 2 at the minimum and the SM minimum. 



where X t = A t — fi cot (3 and m| < m 2 ^ 



are the two eigenvalues of the stop mass matrix 



Mf = 



,2 3 S tu) m Z 



cos 2^ 



m| + m 2 + \s 2 w m 2 z cos 2/3 



(4.20) 



The approximation in (4.19) is valid up to order m^/4m| , but is useful to give a qualitative picture 
of the effect. Nevertheless all the plots and the numbers that are quoted here are obtained with the 
full MSSM one loop result [43]. 

The h/yy vertex correction can be computed rescaling the contribution to the gluon-gluon ampli- 
tude with the appropriate factors of N c , Q 2 and t c (the Dynkin index of the color representation). 
Since the SM h — > 77 loop is dominated by the W boson contribution, the two corrections are opposite 
in sign, with roughly a factor of 4 difference in magnitude 



r G = l + 5r G r 7 w 1 - 0.276r G . (4.21) 

The current data point to \tq\ w 0.9 and |rJ ~ 1.4 as can be seen in the left panel of Figure 5. This 
can be realized only in a small corridor of the stop plane depicted in the left panel of Figure 6, that is 
mostly outside of the 10% tuning region and requires the presence of a very light stop mj < 200 GeV 3 . 
In this corner of parameter space, the stops' contribution to the gluon fusion rate is opposite in sign 
and roughly double in magnitude with respect to the top quark. This region is not only tuned by more 
than 1% if we demand that the gluino be heavier than 1 TeV, but can also generate color breaking 
minima [63]. 

The preferred area with |r 7 1 > 1.1 is rather narrow since the lightest stop mass varies fast for values 



of X* and 



Q3 



that give an r 7 enhancement. This is reflected in a correspondingly rapid 



variation of 5r G as a function of the soft supersymmetry breaking parameters. 



s Allowing for larger values of Xt it is possible to accommodate large h 



stop, for instance we could have At ~ 10 TeV, rrij 



400 GeV and still <5r* 



> 77 enhancements also with an heavier 
i 0.5 [62]. 
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Loop Effects Top Partner 




Figure 5. Right: \ 2 contours at the 68% and 95% confidence levels in the (rc, r 7 ) plane. All other parameters 
are fixed to their SM value. Left: A\ 2 vs. rc for the stops fit, with 95% and 99.73% heights marked in red. 
In gray we show the offset between the \ 2 a t the minimum and the SM minimum. 



The results of the fit are presented in the right panel of Figure 6. The data currently prefer the 
narrow corridor described above. The 3er contours instead lie both in the deeper minimum around the 
large X t region and in the high region where stop effects decouple. This second area corresponds 
to the shallower SM-like minimum of the fit shown in the right panel of Figure 5. 

Overall the results of the fit are in tension with the requirement of 10% tuning, but point to a 
region directly accessible at the LHC. However a good fraction of the preferred masses falls where the 
sensitivity of current searches vanishes: 160 GeV < m t ~ < 220 GeV [64], possibly motivating additional 
efforts. 

4.5 Charginos 

At least one chargino should be light m ± < 300 GeV in view of naturalness considerations (see 
Section 2), possibly giving an observable deviation in the h — > 77 rate. However, the effects can be 
decoupled by raising Mi to 0(1 TeV), without introducing more than a 20% tuning [16]. Nonetheless, 
other than naturalness, there are further motivations to study chargino effects on Higgs couplings. 

It possible that these deviations are the only observable effect of the MSSM at the LHC. Many 
models have the gauginos and higgsinos as the only light particles, with the possible addition of gluinos, 
as is the case of split supersymmetry [65, 66] and models that address the moduli problem and dark 
matter [67-69]. The most studied scenarios typically decouple either gauginos or higgsinos [70, 71], 
which would make the effects on the Higgs couplings vanishingly small, but it is possible to obtain a 
viable dark matter candidate together with unification keeping both species light [72]. Furthermore 
they are a reasonable proxy for a member of a SU(2)l multiplet providing a dark matter candidate. 
The mass splittings in the multiplets are expected to be of the order of a few hundreds of MeV [73] , 
thus leaving open only monojet and monophoton searches at 14 TeV. These analyses may be the only 
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Figure 6. Left: Regions in the stop plane (m^ = m Q 3 = mu 3 ) where the h — > 77 rate is enhanced with 
respect to the SM. Black contours give the lightest stop mass in GeV. In orange percent fine-tuning (A" 1 ) 
contours, as defined in equation 2.2. Right: Shaded areas correspond to 95% (blue) and 99.73% (light purple) 
confidence levels in the (Xt, mj) plane. The black lines are contours of constant mi in GeV. In orange percent 
fine-tuning (A -1 ) contours are shown, as defined in Eq. 2.2. In both panels the gray area corresponds to the 
lightest stop becoming tachyonic. 



viable path to direct exclusion/discovery also in the case of compressed SUSY spectra [74]. Their 
sensitivity is currently limited only to colored particles [75-77] and it was estimated that to probe 
the electroweak production of particles with masses up to 200 GeV, 300 fb -1 at 14 TeV will be 
necessary [78]. Therefore it is worth studying the chargino effects on the /177 coupling in the MSSM, 
that were previously considered also in [72, 79-81]. 

In the MSSM, the chargino contribution is bounded by: —0.3 < 6r* < 0.13. It might be larger 
in the NMSSM [82] or in Z?-term extensions [83], but we do not consider this possibility here. Both 
limits come from direct chargino searches and can easily be derived by scanning (-M2, fi, tan/3) space 
with the LEP constraint rn^± > 94 GeV [84]. The bounds are saturated when tan/3 = 1, restricting 
tan/3 > 2 gives -0.2 < Sr^ < 0.1 [15]. 

It is thus clear that charginos alone can hardly explain the observed h — > 77 excess. The results 
of the fit in Figure 7 roughly point to fiM 2 > sin2/3m^ > where the /177 coupling is enhanced 



± sin 2/3 

7 tiM 2 -m^ sin 2^ 



However the allowed region of the (/x, M2) plane, even for tan/? = 1, is at least 2a away from the 
best fit to the data, while it disappears as tan /3 grows and the effect decouples. Note that the leading 
term in a m^/m 2 ± expansion (equation 4.22) ceases to be a good approximation in most of the region 
preferred by the fit, where one of the charginos can be extremely light 20 — 50 GeV and already ruled 
out by LEP. As in the stops case, all the plots and the numbers that are quoted were obtained with 
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Figure 7. Contours of the \ 2 at 68% and 95% C.L. in the (M2, n) plane for tan/3 = 1 (left) and tan/3 = 5 
(right). The black lines are contours of constant m ± in GeV. In gray is shown the region excluded by LEP 

m ± < 94 GeV. 

"1 



the full MSSM one loop result [43]. 
4.6 Five parameters fit 

In the previous sections we fit the data to simplified scenarios in which all the Higgs rates could be 
expressed in terms of one or two parameters. This gave an idea of the level of agreement between the 
data and some physically motivated corners of the parameter space of natural supersymmetry. In spite 
of the fact that only three channels per experiments have errors below the 50% level (namely 77, ZZ 
and WW untagged) we still find it interesting to explore the more general case in which r 7 , re, r&, r t 
and Ty all play a role. This is a small modification of the four parameter natural MSSM, inspired by 
the possibility of adding a new singlet, and comes closer to approximating a fit with all couplings left 
to float. The only difference with respect to a four parameter fit resides in the fact that r t is virtually 
unconstrained and not artificially limited by its relation with ry . These theory inspired exclusions can 
be applied also outside of the framework of natural supersymmetry and to facilitate possible attempts 
we show profiles of the % 2 for the five couplings in Appendix A. 

In Table 2 we show the corresponding one dimensional 95% C.L. intervals, obtained by treating 
the other parameters as nuisances. The errors on the single couplings vary from 30% to 100%, which 
alone is not enough to lose all hope of constraining the parameter space of natural supersymmetry. 
However the impact of leaving all couplings to float is strong for charginos and the whole plane becomes 
accessible. Similarly for the case of stops, where the contribution in a large part of parameter space 
can be compensated by a shift in r t . 

On the other hand not all sensitivity is lost on tree-level couplings. For instance A35 = A^| SSM 
still gives a preferred region with mg > 250 GeV and A = 1 implies m# > 320 GeV, as it is shown 
in Figure 8. Similarly the tree-level statements about singlet mixing and the decoupling limit are 
meaningful. Repeating the exercise in Section 4.3, but profiling the full five dimensional \ 2 we obtain 
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2HDM 




Figure 8. \ 2 contours corresponding to the 68% and 95% confidence levels in the (A35,m_ff) plane. In the 
case of ri constrained by the five parameters fit (rG,r~/,rv ,rt,rt). The solid purple line corresponds to the 
tree-level value of A35 in the MSSM. The band in yellow covers the possible values of A35 in pure D-term 
models (defined in section 4). The blue dashed line runs through the best fit points. 

I cos 0| < 0.8 at 95% C.L. very close to the result of the purely tree-level fit. This is an indication 
that both the tree-level fit and the five dimensional one are dominated by the lower bound on ry that 
for a ~ /3 — ir/2 4 becomes ry ss sin</>. The bound on cos</> becomes much stronger away from the 
decoupling limit as can be seen in the right panel of Figure 9. 

We can also extract information on the level of decoupling of the Higgs sector as a function of 
tan/?. In Figure 9 the 95% C.L. contour is plotted in the (£, tan/3) plane. The key message is that 
large tan /3 is allowed only for values of £ close to decoupling, where corrections to r& ps 1 — £ tan (3 

4 Note that this is the a defined in Equation 4.17. 



Table 2. Confidence intervals for the five parameters that encode natural SUSY predictions for Higgs rates. 
All down-type couplings scale with r& and alll up- type ones with r t . All other couplings not present in the 
table are fixed to their SM value. 
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Figure 9. Left: The 95% C.L. region in the (£, tan/3) plane, obtained profiling the full five dimensional x 2 - 
The dashed contours indicate r\, for cos 2 <j> = 0. Right: The 95% C.L. region in the (£, cos 2 (j>) plane, obtained 
profiling the full five dimensional \ 2 ■ The dashed contours correspond to different values of rv- 

and r t ~ 1 + £/ tan/3 are small. This is reflected in the (m^jtan/?) exclusion discussed in section 4.2 
that is even competitive with direct searches. Many of these bounds apply to a vast class of complete 
theories and indicate that in the MSSM and many of its motivated extensions, tree-level effects in 
the Higgs sector are already strongly constrained, mainly by the measurements of the h — > WW and 
h — > ZZ rates that in our setting are always below their SM value unless £ = and cos = 0. 

5 Conclusion 

So far most of the effort has been devoted to simplified settings or complete models in which only one 
or two parameters enter Higgs rates. We set upon the more ambitious task of constraining a realistic 
theory using only naturalness to limit the number of parameters entering the game. Some of these 
ideas were treated already in [15, 38], but an application to the data was still missing. 

We fit Higgs measurements, first treating independently the tree-level and loop effects, and then 
performing a five dimensional fit that parameterized complete natural theories such as the NMSSM, 
the MSSM and its D-term extensions. Considering a generic type-II 2HDM we found a bound on the 
heavy CP even Higgs mass run > 370 GeV in theories with an approximate PQ symmetry. We also 
obtained a strong preference for the decoupling limit: |£| < 0.1, and in theories with an extra singlet a 
robust constraint on its mixing. When taking into account only loop effects we showed that the data 
prefer a very light stop (< 200 GeV) with large mixing, while charginos in the MSSM alone can not 
explain the current h —¥ 77 enhancement. 

With the full five parameter fit we found that loop-level statements cease to be valid. However we 
were still able to draw interesting conclusions on tree-level mixings in the Higgs sector. We found that 
the mass bound in a type-II 2HDM with moderate to large tan j3 is relaxed to mji ^ 250 GeV and so 



is the upper bound on singlet mixing, | coscj)\ < 0.8. However there is still a strong correlation between 
the size of tan/3 and the vicinity of the theory to the decoupling limit. This was translated into 
an exclusion in the (rriA, tan /3) plane of the MSSM, that is competitive with direct searches. These 
bounds hold in a large class of complete theories, where loop corrections to the Yukawa couplings are 
small. 

The comparison between the constraints discussed above and LHC searches indicates that there 
is still no tension between the two sets of measurements in natural supersymmetry and that direct 
searches are still the best chance of finding supcrpartners, barring extreme configurations with little 
visible and invisible energy in the event. The only exception are extra states in the MSSM Higgs sector, 
that at moderate tan/3 may be seen first as a deviation in Higgs couplings, rather than produced 
directly at the LHC. 
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A Results of the five dimensional fit 

This appendix collects the one dimensional profiles of the \ 2 obtained from the five dimensional fit 
described in section 4.6. In Figure 10 we show them with the 68% and 95% intervals marked in orange. 
The details of the profiling are discussed in section 3. r t , entering only at loop level in well measured 
rates, is essentially unconstrained due to the compensating effect of tq and r 7 . 
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Figure 10. A\ 2 profiles for the five parameters fit (re r-y, tv, rj,, rt). Heights corresponding to the 68% and 
95% confidence levels are marked in red. 
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